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ABSTRACT 

A simple model for variability in relativistic plasma outflows is studied, in which 
nonthermal electrons are continuously and uniformly injected in the comoving frame 
over a time interval At. The evolution of the electron distribution is assumed to 
be dominated by synchrotron losses, and the energy- and time-dependence of the 
synchrotron and synchrotron self-Compton (SSC) fluxes are calculated for a power-law 
electron injection function with index s = 2. The mean time of a flare or pulse 
measured at photon energy E with respect to the onset of the injection event varies 
as E~^^'^ and E~^^^ for synchrotron and SSC processes, respectively, until the time 
approaches the limiting intrinsic mean time (1 -|- z) At /{I'D), where z is the redshift 
and V is the Doppler factor. This dependence is in accord with recent analyses of 
blazar and GRB emissions, and suggests a method to discriminate between external 
Compton and SSC models of high-energy gamma radiation from blazars and GRBs. 
The quantitative behavior of the X-ray spectral index/flux relation observed from BL 
Lac objects can be explained with this model. This demonstrates that synchrotron 
losses are primarily responsible for the X-ray variability behavior and strengthens a 
new test for beaming from correlated hard X-ray/TeV observations. 



Subject headings: BL Lacertae objects: general — galaxies: jets — gamma rays: bursts 
— radiation mechanisms: nonthermal 



1. Introduction 

Space-based observatories such as ASCA, RXTE, Beppo-SAX and the Compton Gamma 
Ray Observatory, and ground-based air Cherenkov telescopes such as the Whipple Observatory 
and HEGRA, have gathered data from blazars and GRBs of sufficient quality to allow combined 
temporal and spectral analyses. A general trend is becoming apparent in blazar studies (see 
Wagner 1997 and Shrader & Wehrle 1997 for reviews): variability of optical and > 100 MeV 
emission in flat spectrum sources such as 3C 279 (e.g., Hartman et al. 1996) and BL Lacertae 
(Bloom et al. 1997), and of X-ray and TeV emission in the BL Lac objects Mrk 421 and Mrk 501 
(e.g., Macomb et al. 1995; Buckley et al. 1996; Catanese et al. 1997) appears to be temporally 
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correlated, implying that the same population of electrons produces both the optical/X-ray and 

the 100 MeV - TeV emission. In the specific case of the 1994 May flare of Mrk 421, ASCA 
observations (Takahashi et al. 1996) also show that the time lag of photons in the 0.5-2.0 keV 
range relative to 5 keV X-rays varies as £'(keV)~^/^ — 5~^/^. Moreover, the Mrk 421 flare data 
follow a well-defined trajectory in a spectral index/flux display, and this behavior is also found 
in several other BL Lac objects (e.g., OJ 287, Idesawa et al. 1997; PKS 2155-304, Sembay et al. 
1993) 

Recent analyses of GRB data show that the light curves tend to become narrower at higher 
energies (Fishman et al. 1992), and that the decaying phases of GRB light curves are generally 
longer than the rising phases (Link, Epstein, & Priedhorsky 1993). Fenimore et al. (1995) showed 
that the autocorrelation function of 45 bright BATSE GRBs, aligned at their peak fluxes, displays 
an energy dependent width oc E~^-^. In the specific case of GRB 960720, which displays a single 
well-defined pulse, an energy-dependent duration oc E~'^ '^^ is found (Piro et al. 1998). 

The radio/optical continuum of blazars is thought to be nonthermal synchrotron radiation, 
and blast wave models of GRBs also indicate that the hard X-ray/soft gamma-ray continuum 
from GRBs is produced by the same process (e.g., Meszaros, Rees, & Papathanassiou 1994; Tavani 
1996; Waxman 1997; Vietri 1997). The origin of the high energy gamma-ray emission in blazars is 
generally thought to originate from Compton processes, though it remains unclear whether internal 
synchrotron photons or photons produced outside the jet represent the dominant soft photon 
source. The high-energy radiation from GRBs might also originate from Compton-scattering 
processes (see, e.g.. Hurley et al. 1994; Meszaros et al. 1993). The likelihood that synchrotron 
emission from cooling nonthermal electrons can explain some of the previously mentioned trends 
has been discussed in various approximations (e.g., Tashiro et al. 1995; Takahashi et al. 1996; 
Tavani 1996; Kazanas, Titarchuk, & Hua 1998; Kirk, Rieger, & Mastichiadis 1998), but without 
a consideration of the associated SSC emission and a full treatment of the relativistic boost and 
energy normalization for the injected nonthermal electrons. 

In this Letter, we analyze the simplest possible model that retains the essential physics of 
nonthermal electron injection and cooling in relativistic plasma outflows. Even so, it goes a long 
way towards explaining the observed trends, and in addition suggests a fruitful avenue of research 
which could help to discriminate between models of high-energy blazar emission, to determine the 
origin of the high-energy radiation in GRBs, and to strengthen a new beaming test for blazar 
radiation (Catanese et al. 1997). 

2. Analysis 

Assume that nonthermal electrons with Lorentz factors 71 < 7j < 72 are injected uniformly 
throughout a spherical emission region with radius /, and that the electron injection spectrum 
maintains a constant amplitude and power-law form with injection index s between comoving 
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Ne{li,ti) = ,_2-s 2-Ll 0(7i;7i,72) Q{ti;ti,h + At) , (1) 



times h < t < ti + At. Note that causality restricts At to be no shorter than ^ l/c, though 
the injection duration could be much longer if the time scale At is set by some intrinsic process 
which energizes the nonthermal electrons on a longer time scale. If the total energy injected in 
nonthermal electrons during the injection event is £^e(ergs), then the injection function at time t, 
is given by 

(2 - s)£, 
(72"'-7?-^)mec2At 

where the Heaviside function Q{x; a,b) = 1 if a < x < b and Q{x; a,h) =0 otherwise. Equation 
(1) is readily generalized to time-varying injection events and, given some work, to cases where 
the injection occurs nonuniformly throughout the volume of the plasmoid. 

The synchrotron energy-loss rate of isotropic relativistic electrons in a randomly oriented 
constant magnetic field with mean field strength iJ(Gauss) is given by the well-known expression 
—7 = 1^07^, where = ca^H'^ j {^-KnieC^) = 1.29 x 10~^ii"^ s~^. Here we assume that synchrotron 
losses dominate adiabatic losses and Compton losses. The former condition holds for electrons 
with 7 ;^ 8 X 103/?exp/[i^^(G)(V3 x 10^^ cm)], where /?exp is the expansion speed in units of c. The 
latter condition requires that the photon energy density in the comoving frame is less than the 
magnetic field energy density. 

When these conditions hold, the synchroton energy-loss equation is easily solved to give 
7(t) = [7"''^ + VQ{t — ti)]~^ The evolving electron distribution function at time t is therefore given 

by 

i-t /■rmn{t,/\t) 

,(7; t) = / dti Neij, ti) = Xe7"' / dti [7-1 - Mt - ti)Y-^ , (2) 

JO Jmax(OJ-r) 



IS 



where is the coefficient proceeding 7^^* in equation (1) and T = T{'j) = Uq {'y~^ — 72 
the time scale for an electron injected with 72 to reach Lorentz factor 7 through synchrotron 
losses. The expression on the rhs of equation (2), valid for 71 < 7 < 72, is obtained by noting 
that \d'yi/d'y\ = 7^/7^. We have also set ti = without loss of generality. It is straightforward to 
generalize equation (2) for particle escape (see Tashiro et al. 1995). 

We also assume that there is a separation between the acceleration/injection and cooling 
processes. This separation holds when the cooling time scale is much longer than the acceleration 
time scale, but can limit the highest energies to which electrons are accelerated during the main 
portion of a flare, as shown by Kirk ct al. (1998). Particle acceleration can also influence spectral 
and temporal evolution long after the main flaring behavior has ended, though its effect is probably 
negligible if the flares exhibit sharply defined pulses which return to their quiescent levels at high 
energies. 

The spectral flux S(ergs cm~^ e",^^^) observed at energy E = mgC^eobs and observer time 
tobs from a plasmoid at redshift z and luminosity distance cZl is given by 

V^{1 + z) 

'S'(eobsj ^obs) ^obs) = — A-jrd^ — ■ (3) 
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s.n(M) = ^(:^)^/^Ar4(^)V^;.] (4) 



Here fiobs = (/"obsj^obs) specifies the observer direction, where ^obs = arccos /iobs is the angle 
between the jet axis and the direction to the observer, and J(ergs s~^ e~^) is the spectral emissivity 
integrated over the volume of the plasmoid. Equation (3) is valid for radiation which is emitted 
isotropically in the comoving frame, and applies to synchrotron and SSC processes (see Reynolds 
1991 and Dermer 1995 for cases where this does not hold). The expressions e = (1 + z)eobs/^ 
and dt = Vdtohs/iX + ^) relate the comoving and observed photon energies and the comoving and 
measured differential time elements, respectively, where the Doppler factor T) = [T{1 — S/iobs)]~^ 
and B = {1 — l/r^)-*^/^, and T is the Lorentz factor of the plasmoid. We assume that F remains 
constant throughout the flaring event, which implies that if if swept-up matter encrigizes the 
plasma, then the mass of the swept-up matter is -C A^/r, where Ai is the total mass of the 
particles in the plasmoid. A self-consistent treatment of plasmoid dynamics when the inertia of 
external matter swept up by the plasmoid cannot be neglected has been treated by Chiang & 
Dermer (1998). 

In the (5-function approximation for synchrotron emission, which is a good approximation 
away from the cutoffs of the synchrotron spectrum produced by the endpoints of the electron 
distribution, we have 

2caTUH 

(e.g., Dermer, Sturner, &; Schlickeiser 1997). The term uh = H'^/Stt is the magnetic-field energy 
density, and en = H/4A14 x lO^^G is the dimensionless electron plasma frequency. 

After substituting equation (4) into equation (3) and using equation (2) for the electron 
spectrum, we find for the case s = 2 that 

5'syn(eobs; ^^obs, ^obs) = -fTT {min(t. At) - max[0, t - TU — )]} , (5) 

^obs 

where the coefficient Cgyn = P^/^(l + z)^^"^ KgCaTUH / {QTrd'j^e]^'^). The case s = 2 is arguably 
the most interesting case since, according to simple shock acceleration theory, first-order Fermi 

acceleration by a strong shock produces a particle spectrum with s = 2 in nonrelativistic 
monatomic gases, provided that nonlinear feedback of the accelerated particles can be neglected. 
More general injection functions can be considered in more detailed treatments. 

The volume-integrated SSC emissivity in the ^-function approximation is given by 

Jssc(6,i) = ^^5^^ / de'e'-iiVe(J-;t)iVe(J4;t) (6) 

(Dermer et al. 1997). The upper limit on the integral restricts the scattering to the Thomson 
regime. After substituting equation (6) into equation (3) and making use of equation (2) for the 
electron distribution, one obtains, again for the case s = 2, the result 

Cssc ,_ , fV 

'S'ssc(eobs;f^obs,*obs) = — / dee {min(t. At) - max[0, i - r( J — )]} 



-5- 



X {min(t, At) - max[0, t - T{^)]} . (7) 



The coefficient Cssc = CsynKeCfT/ {2ttI'^). 



3. Results and Discussion 

Figure 1 shows the time-dependence of the synchrotron (Fig. la) and SSC (Fig. lb) spectral 
fluxes for a variety of different observing energies, obtained by numerically solving equations (5) 
and (7), respectively. Because S oc e^^l"^ from uncooled electrons injected with s = 2, the spectral 

1/2 

flux is multiplied by ej^^ for clarity of presentation on a linear scale. In this calculation, we let 

= 0.1 Gauss, 7i = 10, 72 = 10^, z = 0.1, T> = 10, and At = 10^ s, as might be appropriate for a 
blazar flare. 

We wish to determine the mean time (t(eobs)) of the radiation observed at different energies, 
measured with respect to the onset of the injection event. This is given by the expression 

/— 00 '^^obs ■ ^obs ■ 'S'(£obs) ^obsj^obs) 
I— 00 ^^obs ' '^'(Cobs) ^obs) ^obs) 

Note that as defined here, a Heavisidc (or boxcar) function has a mean time equal to one-half its 
temporal duration. Equation (8) can be easily generalized for for higher moments of the temporal 
proflle giving, for example, the FWHM duration of a flare. 

After substituting equation (5) into equation (8), one obtains two cases depending on whether 
T < At of T > At. When 72 — > 00, the same result is found in both cases, namely 



=^iA«+-„-y (^1^1. (9) 

Equation (9) provides a convenient expression for fitting energy-dependent time-lag data, such 
as the Mrk 421 flare (Takahashi ct al. 1996), or the energy-dependent GRB widths measured by 
Fenimore et al. (1995) and Piro et al. (1998), though the appropriate moment analysis should be 
used in more detailed treatments. Equation (9) indicates that the energy-dependence of the mean 
time varies more slowly than E~^-^ when the light-crossing time l/c or duration of the energization 
event are comparable to or longer than the comoving electron cooling time scale. Fitting high 
quality data to equation (9) to determine the photon energy eobs where the two branches of the 
expression intersect, and using the shortest variability time scale St^f = (1 -|- z)At/'D observed at 
Cobs ^ fobs; "we derive the plasmoid magnetic field 

^'°>^°-«t ^E(4,;i(,^)P >-/^ (10) 

Causality implies that I ^ cD5f^^ /{I + z). It will be important in future studies to understand 
how equations (9) and (10) are modified for a non-spherical plasmoid geometry and (especially for 
GRBs; see, e.g., Fenimore, Madras, & Nayakshin 1996) a blast wave geometry. 
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Expressions similar to equations (9) and (10) have been noted (e.g., Tashiro et al. 1995; 
Takahashi et al. 1996; Tavani 1996; Buckley et al. 1998; Catanese et al. 1997; Kazanas et al. 1998), 
but this treatment provides a precise fitting function for analysis of data when s = 2 and yields 
a generalization for arbitrary values of s. In the paper by Catanese et al. (1997), moreover, it is 
noted that correlated X-ray/TeV data imply an upper limit on H because the electrons producing 
the highest energy synchrotron emission have Lorentz factors 2 x 10^£^c(TeV)(l + z)/V, where 
-Bc(TeV) is the measured energy in TcV of the highest-energy gamma-rays. Synchrotron emission 
correlated with the TeV flux requires that the electrons radiate in a magnetic field at least as great 
as H{G) = lleobs,syn^/[-Ec(TeV)(l + z)], where 

^obsjsyn is the measured dimcnsionless energy 
of the highest energy synchrotron photons produced by the electrons which produce the TeV 
radiation. When compared with the value of H inferred through equations (9) and (10), we obtain 
an expression for the Doppler factor, given by 

^£(keV)(«X")''^(«ob.,^)=''' 

A lower limit to D is obtained if the TeV flux does not exhibit a clear cutoff due to the high-energy 
cutoff in the the electron distribution function. 

The inset to Fig. la shows the energy dependence of the mean time {t{E)) for the synchrotron 

and SSC processes. As can be seen by examining Fig. lb, the SSC emission decays more slowly 
than the synchrotron process, yielding an energy dependence oc E~'^^^ until the shortest mean 
time (1 + z)At/2'D is reached, which occurs at much higher energies than for the synchrotron 
process. Thomson scattering of external photons, having an energy loss rate of the same form 
as the synchrotron energy loss rate, also gives (t(eobs)) cx^ E~^l'^ for s = 2. When more sensitive 
blazar 7-ray observations become available with the upcoming INTEGRAL and GLAST missions, 
the energy-dependence of the mean time or duration of blazar flares can be used to determine 
whether SSC or ECS processes dominate in specific sources or between different source classes. 
This test might also be possible for bright TeV flares from BL Lac objects given the steadily 
improving sensitivity of air Cherenkov telescopes. More detailed studies need to be performed 
which take into account Klein-Nishina effects and general values of s. 

Figure 2 shows the evolution of the synchrotron spectrum at different observing times, 
using the same parameters as in Figure 1. At a fixed photon energy, the flux rises due to the 
accumulation of injected nonthermal electrons until cooling plays an important role in depleting 
the electron spectrum (cf. Dermer & Chiang 1998). The effects of cooling are seen earliest at the 
highest photon energies, and cause a softening of the spectrum. Consequently, the spectral flux at 
a given photon energy displays a hard spectrum while its intensity is increasing, since the electrons 
that are radiating this emission have not yet felt the effects of cooling. When synchrotron cooling 
becomes important, the flux begins to fall rapidly. This behavior is shown by the solid curve (1) 
in the inset to Figure 2, which displays the clockwise evolution of the 2-10 kcV spectral index as 
a function of a quantity proportional to the 2 keV spectral flux. In reality, there will be some 



-7- 



level of background which will be reached. This is crudely modeled here by adding an underlying 
spectral component with spectral flux 5bg(eobs) = -^bgCobs^' '^here K\^g = 1 and 5 x 10'^ ergs cm~^ 
for the dashed curve (2) and dotted curve (3), respectively. 

The behavior illustrated in the inset to Fig. 2 is in qualitative agreement with evolutionary 
tracks of the spectral index/flux observed from some blazars, as noted in the Introduction (see 
also Kirk et al. 1998). Precise fitting of such tracks will require extending this model for general 
injection indices and for various forms for the background emission. The well-known hard-to-soft 
evolution of GRB pulses (e.g., Norris et al. 1986) could be a manifestation of this effect, and the 
approach outlined here can also be used to analyze the fluence dependence of the peak of the GRB 
uFy spectrum (Liang & Kargatis 1996). 

In summary, a very simple model has been presented for the observed nonthermal synchrotron 
and SSC emission emitted by cooling nonthermal electrons which are injected over a comoving time 
interval Ai into plasma with relativistic bulk motion. In spite of the model's simplicity, empirical 
trends that have become better defined through recent combined temporal and spectral analyses 
of data from blazars and GRBs are qualitatively understood. Straightforward generalizations, 
necessary for detailed fits to data, were indicated thoughout the analysis and will be treated in 
future studies. 

I thank Ed Fenimore and Michael Catanese for conversations and questions that focused 
my attention on this problem. Comments by J. Chiang and the referee are acknowledged. This 
research was supported by the Office of Naval Research and the Compton Gamma Ray Observatory 
Guest Investigator Program. 
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Fig. 1. — Synchrotron (Fig. la) and SSC (Fig. lb) fluxes as a function of time at different photon 
energies. Nonthermal electrons with injection index s = 2 are injected for At = 10^ seconds into 
a spherical plasmoid at redshift z = 0.1, and with Doppler factor V = 10 with respect to the 
observer. Other parameters are given in the text. The mean time of the fluxes as a function of 
observer energy, defined by equation (8), are shown in the inset. 
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Fig. 2. — Synchrotron fluxes as a function of observing energy at different observing times 
iobs(s) = [(1 + z)At/'D)] X 10^~'^^/^, using the same parameters as in Figure 1. The curves are 
labeled by j. The inset shows the tracks followed by the 2-10 keV spectral index versus a quantity 
proportional to the 2 keV spectral flux. Various levels of background fluxes, as described in the 
text, are included in the calculations of the dashed and dotted curves. 



